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1. Introduction 


In recent years pioneering contributions have been made in the field of rigorous 
computer-assisted results for dynamics of dissipative PDEs [21151151151151151 nirlfTTinUl 
(201121]. The methods exploit the smoothing property of the system to apply either 
topological or functional-analytic tools. However, little attention has been paid 
to apply these methods to other types of evolution PDEs, such as the ones with 
tail of saddle type. In such problems we need to deal with an infinite number 
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of strongly repelling and strongly attracting directions. From the point of view 
of topological methods this situation is just as good as the dissipative one - for 
example Theorem 2.3 in m is formulated in a way that is readily applicable to 
finding equilibria of these systems. Our goal in this paper is to take this approach 
one step further and use a topological tool of periodic isolating segments to prove 
the existence of periodic in time solutions in a nonautonomously perturbed equation 
of such type. 

Our example will be the Boussinesq equation [4] but without much difficulty the 
methods can be applied to produce similar results in other systems with an indefinite 
tail of saddle type. 


1.1. The forced Boussinesq equation 


We consider the following second order nonlinear equation perturbed by a time- 
dependent forcing term: 


'Utt — 'Uxx T (3'U‘xxxx + ) xx ~t“ . 


( 1 ) 


On u and / we impose periodic and even boundary conditions and a zero-average 
condition in x : 


u(t,x + 2ir) = u(t,x), (2) 

u(t,—x) = u(t,x), (3) 

p2n 

/ u(t , x)dx = 0 , (4) 

Jo 

f(t,x+ 2 tt) = f(t,x), (5) 

= f(t,x), (6) 

p2ir 

/ f(t,x)dx = 0. (7) 

Jo 

For fi > 0 the unperturbed equation 

^tt — Uxx T (J'Uxxxx + fr(u ) xx (8) 


is the “bad” Boussinesq equation and was derived by Boussinesq [4] as a model for 
shallow water waves. The equation is famous for its ill-posedness. Indeed, when 
looking at its linear part 

Utt. = u xx -f (3u xxxx (9) 

one can observe a rapid growth in high Fourier modes for almost all initial data, hence 
a consequent loss of regularity of the solution. This is a significant complication in 
the numerical analysis of ([8|l. since slightest perturbations of the initial problem 
can produce a totally different behaviour at output. Because of that, regularized 
versions of the equation were considered in numerical studies m- Solutions to 
the equation (JHJ) were also obtained analytically [TQj and by the inverse scattering 
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Tab 1.: Bounds on e and on the norms of periodic solutions and their time deriva¬ 
tives for / £ J 


p 

6 

IKV)II L 2 < 

IKMIloo < 

IMV)|| L 2 < 

IMv)|| c ° < 

1.5 

[-0.05,0.05] 

0.28862115 

0.12440683 

0.24610779 

0.2143523 

1.75 

[-0.1,0.1] 

0.41504192 

0.1820825 

0.39340084 

0.42461205 

2.5 

[—0.3,0.3] 

0.84724825 

0.38676747 

0.96839709 

1.4795576 


Tab 2.: Bounds on e and on the norms of periodic solutions and their time deriva¬ 
tives for / € Jq? 


P 

6 

\\ U (t,-)\\ L 2 < 

IKt,-)llco < 

IMV)|| L 2 < 

IMv)|| c ° < 

1.5 

[-0.05,0.05] 

0.29831987 

0.13194161 

0.25703095 

0.24099758 

1.75 

[-0.1,0.1] 

0.43198386 

0.19524766 

0.41478653 

0.47720834 

2.5 

[-0.3,0.3] 

0.88825406 

0.41784158 

1.0309512 

1.637095 


method m ■ Our approach is different; we analyze the direction of the vector field 
on certain subsets of the phase space, and by a topological method we deduce the 
existence of smooth, periodic solutions. 

Here is an example result illustrating our method. We define families of functions 

■?>={/ : f(t, x ) = 2/i(f)cosa:, 

/i continuous and r-periodic, |/i(t)| < 1 Vi} , 

R (10) 
Ft = U'■ f(t,x) = 2^2 f k (t) coskx, 
k =1 

fk continuous and r-periodic, \fk(t)\ < 1 Vfc,<} . 

Then the following theorem holds. 

Theorem 1. For a = 3, for all t > 0 and all f(t, x ) £ U Ff? , and for values of 
P and e given in Tables Q] and[H there exists a classical r-periodic in time solution 
to CD, subject to conditions 0, 0 and 0- The solution exists in the vicinity of 
0 and the bounds on its L 1 2 and C° norms and the norms of its time derivative are 
given in Tables 171 The solution and its time derivative are C 4 and C 2 smooth 

in x, respectively. 

Observe, that 0 is a constant in time solution of the unperturbed system 0, 
hence the requested r-periodic solution for e = 0. Nevertheless, the method is not 
perturbative. We consider a perturbation problem only because it gives a convenient 
approximation of the periodic solution for |e| ^ 0 small. 

The proof is computer-assisted, that means certain inequalities contained in it 
are verified rigorously by a computer program in interval arithmetics. The program 
source code is available at |J]. From Table [3] and equation (l26l) one can also extract 

1 By bounds on L 2 and C° norms of a function u = u(t,x) we mean upper bounds on 

suPtem IMb 01 La and sup teK ||u(i, 0M c o> respectively. 
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the exact bounds on the Fourier coefficients of the solutions, which we do not give 
here. By following the steps of the proof it will become clear that we can easily 
produce results of the same type for: 

• any parameters cr G M and (3 > 1, 

• any given smoothness s > 5, 

• any forcing of the form f(t , x) = fk(t) cos kx where each fk is continuous 

and r-periodic. 

The periodic solution will exist for e € [— £q, £o], £q small enough, and we can attempt 
to verify an explicit range and obtain a bound for the norm with help of the program. 

Let us finish this section with listing some generalizations. Most of them can be 
adapted from the previous treatment of dissipative PDEs sum] and several would 
require only little effort to be introduced in this paper. However, we consider this 
exposition as a preview of the method. We tried to focus on the key matter, which 
is how to deal with the linear instability of the high Fourier modes in a simplest 
scenario. 

1. We could take /3 from the range (0,1]. Then, the linearized equation [9] possesses 
purely imaginary eigenvalues in its low modes. This would involve conducting a 
finite-dimensional analysis of the higher order terms of the low modes. 

2. We could allow non-zero averages, non-even functions or periodic solutions ob¬ 
tained in the proximity of non-zero equilibria of the unperturbed system. 

3. For the forcing term it would have been enough to assume a sufficiently fast decay 
in the high Fourier terms. 

4. We could consider non-periodic forcings and attempt to prove existence of (not 
necessarily periodic) solutions that exist for all t G R, by techniques from [6]. 

The apparently difficult problems are 

• proving the existence of periodic solutions which are not obtained as pertur¬ 
bations of stationary points, 

• proving dynamics more complicated than a periodic orbit (e.g. chaotic dynam¬ 
ics), 

• proving the existence of periodic orbits in autonomous ill-posed systems. 

We think that to efficiently treat these cases within the framework of self-consistent 
bounds we would need a rigorous integration procedure, akin to the rigorous inte¬ 
gration of dissipative PDEs [3[T9l|20]. Obviously, the ill-posedness is a significant 
issue and it seems that the integration should be combined with an automatic seg¬ 
ment placement in the expanding coordinates. We are currently looking into the 
feasability of this approach. 

This paper is organized as follows. In Section [2] we invoke the general method of 
self-consistent bounds and a result which states that a sequence of solutions for the 
Galerkin projections converges to a solution of the PDE. In Section [3] we present a 
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result of Srzednicki m stating the conditions under which non-autonomous time- 
periodic ordinary differential equations - in our case the Galerkin projections - have 
periodic solutions. In Section 0] we apply these tools to the Boussinesq equation JTJ 
and prove Theorem [I] 


2. The general method of self-consistent bounds 


In this section we recall the general method of self-consistent bounds, as introduced 
in the series of papers 19 5D.I2I • We follow the exposition given in 6j for time- 
dependent systems. 

Let J C K be a (possibly unbounded) interval. We consider a nonautonomous 
evolution equation on a real Hilbert space H ( L 2 in our case) of the form 

d i = p M. (ii) 

We assume that the set of x such that F(t, x) is defined for every t £ J is dense 
in H and we denote it by H. By a solution of (fill) we mean a function a : J' —> H, 
such that J' is a subinterval of J, u is differentiable and m is satisfied for all t £ J'. 
Let I C h d and let H k C H be a sequence of subspaces with dim H k < d\ < oo, 

such that _ 

H = 0 H k (12) 

fee/ 

and Hk s are pairwise orthogonal. We will denote the orthogonal projection onto 
Hk by Ak and write 

a fc := A k a . (13) 

From m it follows that a = a k - 

From now on we will fix some (arbitrary) norm | • | on Z d . For n > 0 we set 

X n ~ 0 H k , 

kel,\k\<n (14) 

W := Xn ■ 


We will denote the orthogonal projections onto X n and Y n by P n : H —> X n and 
Qn'-H Y n . 

Definition 2. Let J Cl be an interval. We say that F : J x H D domF —> H is 
admissible, if the following conditions hold for each i £ lA such that dim Xi > 0: 

• J x Xi C domF, 

• (Pi ° F)\ JxXi : J x Xi —> Xi is a C 1 function. 
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Definition 3. Let F : J x H —» H be admissible. The ordinary differential equation 

^ = P n F(t,p), pel„, (15) 

will be called the n-th Galerkin projection of CD- 

Definition 4. Assume F : J x H —► Ft is an admissible function. Let m, M £ N 
with m < M. A compact set consisting of W C X m and a sequence of compact 
sets {Bk}kGi.\k\>m such that Bk C Hk form self-consistent bounds if the following 
conditions are satisfied: 

Cl For |fc| > M, k € I it holds that 0 £ Bk, 

C2 Let |fc| > to and hk '■= max ae s fe ||a||- Then '}2\k\>m kei < 00 ■ I n particular 
we have 

IE© B k cH. (16) 

|fc|>m 

C3 The function (t,u) —>■ F(t,u) is continuous on J x W © ri|fc|>m-®fc C K. x 
H. Moreover, if we define f k := su P(t,u)eJxW®U keIAkl>rn B k \\ A kF{t,u)\\, then 

^-^|fc|>ra,fc£ J fk < '' 

Given self-consistent bounds formed by W and { Bk}kei.\k\>m > by T (the tail) we 
will denote the set 

T:= J] B k CY m . (17) 

k£l ,\k\>m 

The following theorem is a straightforward adaptation of Lemma 5 from [TS] (see 
also Section 4 in m) to a nonautonomous setting. 

Theorem 5. Let W and {Bk\k^i,\k\>m form self-consistent bounds and let 
{nfejfegN C N be a sequence, such that lim^oo nk = oo. Assume that for all k > 0 
there exists a solution Xk '■ [t \, £ 2 ] —> W © T of 

= p n k p(t)£X nk . (18) 

Then there exists a convergent subsequence lim^ooPfc, = p*, where p* : [ti, < 2 ] —>• 
W © T is a solution of (HU)- Moreover, the convergence is uniform with respect to 
t on [t\, f 2 ] • 


It turns out that it is fairly simple to find self-consistent bounds. In the treatment 
of evolution PDEs such as Kuramoto-Sivashinsky or Navier-Stokes it is enough to 
take tails of the form Bk = {a £ Hk : ||a|| < C/ |fc| s } for s large enough. This will 
also be the case in our study of the Boussinesq equation. 
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3. Periodic isolating segments 


The purpose of this section is to recall a result of Srzednicki [12] on the existence of 
periodic orbits in non-autonomous time-periodic ODEs. For the Boussinesq equation 
this theorem will be used to treat each of the Galerkin projections of the system. 
We consider an ODE 

x = g(t,x), (19) 

where g:Ix R n — > R” is of class C 1 and r-periodic in t. Let L 1 : I" —> R, * = 
1,..., k be C l functions and let r € {0,..., k} be fixed. We define sets 

So ~ {leK": L\x) < 0 Vi = 1,.. 

Sq :={*6R":3i6l.r: L\x) = 0}, 


and put S := [0, t] x Sq, S := [0, r] x 5 0 . 

Definition 6. We call S a periodic isolating segment over [0, r] and S~ its exit 
set iff the following conditions hold 

(51) g(t, x) ■ 'S7L l (x ) > 0 for t G [0, t\, i G {1,..., r} and x G So ■ L z (x) = 0. 

(52) g(t, x) • VL*(x) < 0 for t G [0, t\, i G {r + 1,..., k} and x G So : L’ , (x) = 0. 

Theorem 7 (Theorem 2 in |12j). Let S be a periodic isolating segment over [0,r] 
for g defined as above. If S and S~ are compact absolute neighborhood retracts and 
the difference of their Euler characteristics x(S'o) — x{Sq) is non-zero, then there 
exists a point Xo G int So such that the solution x(t) of (1T51) satisfying x(0) = xo is 
r-periodic in t and x(t) G int So for all t G R. 

In what is below we denote the i-th coordinate of a vector x G R n by Xi. 

Corollary 8. Let So := [a\, 6i] x • • • x [a n , b n \ and g = (gi,..., g n ) be defined as 

above. Suppose that for each i it holds that 

gi(t,x)gi(t,y) < 0 (21) 

for all x,y G So ■ Xi = ai, yi =bi and for all t G [0, t\ . Then there exists a r-periodic 
in t solution of (ED with values contained in int Sq ■ 

Proof. After rearranging the coordinates we can assume that 

g t {t,x)> 0, for x G So : x t G {bi,..., b r , a r +i,. ■ ■, a n } , ^ 

g t (t,x) < 0, for x G S 0 : x t G (ai,..., a r , b r+ i,..., b n } . 

for all t G [0, r] and some r G {1,..., n}. 

Let L l (x) := (xi — ai){xi — bf), i = 1,... ,n. Then So is given by the sequence 
{L 1 } and 

VL l {x) = 2 Xi — ai — bi. (23) 

2 The original paper m uses the notion of ( p, g)-blocks, periodic isolating segments are intro¬ 
duced in later studies Ena as a more general tool. 




Now conditions (SI) and (S2) immediately follow from E^l) . hence S = [0,r] x So 
forms a periodic isolating segment. Since Sq consists of r opposite faces, we have 
X'(S'o) — x(5o") = (—l) r and the assertion of Theorem [T] follows. □ 

Later on we will refer to condition m as the isolation conditions or the isolation 
inequalities. 


4. Application to the Boussinesq PDE 


4.1. The Boussinesq equation in the Fourier basis 


Our goal in this section is to express our PDE in coordinates suitable for application 
of Corollary [8] to the Galerkin projections. For that purpose we express the problem 
in the Fourier basis and diagonalize its linear part. 

By formally substituting u(t, x) = X^kez u k(^) elkx into the Boussinesq equa¬ 
tion 0 we obtain an infinite ladder of second order equations 

u k = k 2 {/3k 2 - l)u k - ah 2 ^ u k iU k _ kl + f k (t), kGZ. (24) 

fciez 

Since u is real and even in x , we have u k = U- k and u k € R for all k € Z. Moreover, 
from 0 we have uq = 0. After these substitutions and rewriting the system as a 
first order system we obtain the following equations 


— Vk •> 


v k = k 2 (/3k 2 - l)uk ~ 2ak 2 ^ u^+kU^ ~ crk 

fci>i 


k -1 

2 E 

k\=l 


u kl u k - kl + ef k (t ), k € N + . 

(25) 

As one can see, the linear part of (E5l) is already in a block-diagonal form. All 
we need is to diagonalize each of the blocks. From now on we assume that /3 > 1. 
After a simple calculation we see that the eigenvalues of the linear part of (1^51) are 

r _ -1 T 

±\/k 2 ((3k 2 — 1) with eigenvectors 1, ±^'k 2 (fdk 2 — 1) , respectively. We introduce 

the variables u j!" and such that 


We have 


Uk 


1 

1 

[ < I 

v k 


V' k 2 {pk 2 - 1) 

-^k 2 {/3k 2 - 1) 

. u k . 


1 

2 

1 

2 


_l_ 

2^/fc 2 (/3fc 2 -l) 

2^/k 2 (/3k 2 -l) 


Uk 

V k 


(26) 


( 27 ) 
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and our equations become 


\Jk 2 (f3k 2 - l)ul 


-^k 2 ((3k 2 - 1 )u k 


ak 2 N k (u) + ef k (t) 
2\/k 2 (/3k 2 - 1) ’ 
_ <jk 2 N k (u) + ef k (t) 
2^/k 2 {/3k 2 - 1) 


(28) 


where 

k -1 

7V fc (u) := -2 ^ ^ u kl u k - kl . (29) 

/Ci>l /Cl — 1 

We remark that it is unprofitable to rewrite the convolutions in the new variables 
as we will eventually estimate these terms. 

The n-th Galerkin projection of (E51) is given by 


u k = \/k 2 (l3k 2 - 1 )u^ + 


u k = -^k 2 {Pk 2 - 1 )u k 


<yk 2 N k ^ n (u) + ef k (t) 

2 ^k 2 {/3k 2 - 1) ’ 

_ ak 2 N ktn (u ) + ef k (t) 

2 ^k 2 {/3k 2 - 1) 


k = 1,..., n , 


(30) 


where 


n—k k —1 

* = 2 ^ ^ Uk 1 +kUk 1 ^ ^ . 

fci>l k±=l 


(31) 


Our next step is to construct a sequence of isolating segments S n for the Galerkin 
projections (l30l) . 


4.2. Construction of periodic isolating segments 


We look for periodic isolating segments S n of the form S" = [0 ,t]xSJ, where 

M n 

S 0 ” = IIK’^] 2 ® II l~C/k s ,C/k s } 2 , (32) 

fe=1 k=M +1 

i.e. the set of n-tuples of pairs (u + ,u~) = {{u k , u k )} k=1 such that u k ,u k £ [u l k ,u k ] 
for k £ 1 and \u k \ < C/k s . For now it is enough to take C £ R + and 

s £ {2,3,...}, however later on we will assume that s is at least 6, to comply with 
condition C3 from the definition of self-consistent bounds. 

Observe that we would like to choose the values of C and s, as well as the first 
M intervals the same for each projection. Therefore we can say that our segments 
are a projection of an “infinite-dimensional segment” given by 

M oo 

n b c/k s ,c/k s } 2 . 

k =1 k—M-\-l 


( 33 ) 
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The elements of Sfi° are sequences of pairs (u + ,u~) = {{u k ,u k )} k L 1 . We denote 
them by the same symbols as elements of Sq but we will always make it clear to 
element of which set we are referring to. 

We would like to choose u l k , u r k , C and s such that the linear part of (l30l) domi¬ 
nates the nonlinear terms and the isolation conditions m hold - at least for suffi¬ 
ciently high modes, for n large enough. The inequalities for the low modes we will 
treat one-by-one with aid of rigorous numerics. 

We assume the bounds for u k to be the same as the ones for u k . As we will see 
later, due to the symmetry of the equations (l30l) . the isolation conditions for both 
u k and u k are given by the same inequalities. 

From estimates in it follows that, for a set Sjf of the form given above, 
there exists a constant D £ M + such that 

sup {|iV fc ,„(ii)| : n> k > M, (u + , u~) £ Sq } < . (34) 

The value of D can be given by an explicit formula, but we postpone its evaluation 
to Subsection 14.31 

Lemma 9. Assume that for some M £ N + , we have fk = 0, k > M and 


c ~~ ° D 1 

M + 1 2(/3 — (Af + l)- 2 ) ' 

(35) 

Consider the Galerkin projection (1301) for n > M and let Sq be given 
for M < k < n and ( u + ,u~) £ Sq the following inequalities hold 

by (l32l) . Then, 

ut> 0 if u t = C/k s , 

(36) 

ii+<0 if u+ = -C/k s , 

(37) 

ii k < 0 if u k = C/k s , 

(38) 

u k > 0 if u k = —C/k s . 

(39) 

Proof. We will prove (l36l) and (l38l). The proof of (l37l) and (l39l) follows by reversing 
the inequality signs. We want 


(40) 

Since uf = C/k s , the above is equivalent to 


C aN kjn (u) ^ 

k s 2 (/3k 2 - 1) ' 

(41) 

By the estimate (1M1) it is enough that 


aDk a D 

> 2(/3fc 2 - 1) “ 2fc(/3 — k~ 2 ) 

(42) 

and the right-hand side is at most • 2 (p~(M+ 1 )- 2 ) • 

□ 

3 We note that the estimates from [20] improve the bound on d34t) to D/k s for 
one we use here is fine enough for our applications. 

some D, but the 
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4.3. Estimates for the nonlinear terms 


In this subsection we provide bounds for the nonlinear terms and compute D. In 
fact we will look for an estimate 


N k , 


< 


D 


k > M 


for 


Nk .max ■•= sup{|7Vfc(«)| : (u + ,u ) e 5“} . 
as it is an upper bound on the left-hand side of (1551) for all n, k : n > k. 


(43) 

(44) 


Let Sfi° be of the form as in (1551) and {u + ,u ) £ Sq°. Recall that u k = u k + u k , 


hence 


Uk £ [2u l k ,2u k ] for k < M , 

2 c 

\uk\ < -— for k > M . 
k s 


(45) 


The nonlinearity consists of terms given by an infinite sum and a finite sum 

N k (u) = —2IS(k) — FS(k) , (46) 

where 

IS(k) := ^ u kl +k'Uk i , (47) 


fc x >i 
k -1 


FS(k) := ^2 u kl u k - kl ■ 


(48) 


k i=l 


These terms, arising from the nonlinearity in the Kuramoto-Sivasinsky equation, 
were estimated in pH] (cf. also [20], Section 8). Throughout the rest of this 
subsection we will denote by u k the whole interval [2u k ,2u k ] and put |ztfc| := 
2max{|it[.|, \u r k \\. 

Lemma 10 (Lemma 3.1 in [H]). For k E {1, ..., M} we have 

M-k M 


IS{k ) C ^2 Wfci+fcUfei +2C' ^2 

k±=l k 1 =M-k-\-l 

4 C 2 


|«fc| 


(k + k i)‘ 




(49) 


+ (k + M + l) s (s-l)M s -^ 1 ’ 1 ^' 
Lemma 11 (Lemma 3.5 in [21]). For k > 2M we have 


FS W C T^T 


2 C ( 2 S+1 


M 


2M+ 1 


5 Z ki 


c 2 2s+1 


k i = l 


C2 S+1 


(2M + l) s+1 (s — 1)M S 


[-1,1] - (50) 


Lemma 12 (Lemma 3.6 in J5Tj). For k > M we have 


2C ( 2C M \ 

lsm c t -M« +1) ( (M + !)-(, -1) + g I- 1 ' ‘1 • < 51 > 
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Following [2D] we give Di , D 2 such that 

\FS(k)\<^, | IS(k)\<^, k > M (52) 

and then set D := Di + 2D 2 . Using Lemmas 1111 and [T21 we get the following formulas 
for Di, D 2 : 


Di(k < 2 M) 

D\(k > 2 M) 

Di 

D 2 


max{fc s 1 |i r S'(fc)|, M <k < 2 M} , 


2 C 


2 s +i 

2M + 1 


M 


m 


C2 2s+1 


k i=l 


C2 S+1 


(2M + l) s+1 (s — 1)M S 


max(Hi(fc < 2M), > 2M)), 

m \ 

) ■ 


2C 


2 C 


M + l \ (M + l) s_1 (s — 1) 


k =1 


(53) 

(54) 

(55) 

(56) 


4.4. Low mode isolation and a procedure for refining the bounds 


We will now discuss the low mode isolation inequalities. Assume, that we found 
MeN and segments Sq , n > M such that the assumptions of Lemma [9] hold. To 
apply Corollary [8] to all Galerkin projections it is now enough to check 


K 

> 0 

if 

+ r 

K = u k . 

(57) 


< 0 

if 

Wfc = 4, 

(58) 

K 

< 0 

if 

Wfc = u r k , 

(59) 

K 

> 0 

if 

s 

»■ 

II 

s 

?r'~~ 

(60) 


for (u , u + ) S Sq and n > M. It is enough to verify 


ak 2 N k {u) + efk(t) 
2 k 2 {/3k 2 - 1 ) 
ak 2 N k (u ) + efk(t) 
2 fc 2 (/3fc 2 - 1) 


(61) 


for all u £ S£°, t £ [0, r] and k = 1,..., M. Recall that Nk(u) = —2 IS(k) — FS(k). 
The term IS(k) is bounded by use of Lemma 1 101 while FS(k) are finite sums which 
can be for example rigorously enclosed by use of interval arithmetics. Therefore we 
can compute an explicit bound 


aN k (u)c[N l k ,N r k ], u£SZT 


(62) 


for each k = 1,..., M. Assume that we are also have some bounds 

ef k (t)c[f l k ,r k \, tel. 


( 63 ) 
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Given these enclosures, inequalities m can be checked easily, in our case on a 
computer using interval arithmetics. Note that there is no guarantee that for given 
e / 0 and given bounds S the inequalities (RJTl) will be satisfied. However, for |e| 
small enough “good” bounds should exist. Since we cannot expect to choose the 
correct values for C and s at first try, we will use an algorithm from m 

(Section 3.3) for refining an initial guess for the bounds. 

Our goal is both to increase s and correct our guesses for the bounds on coordi¬ 
nates where the isolation inequalities do not hold. We iteratively adjust the pairs 
(u k ,u k ), and later C and s, so at each step our new guess is at worst case an equality 
in the isolation conditions. This way the bounds are tight on each coordinate, so the 
nonlinear terms do not contribute much error. Note that the procedure is heuristic 
and we do not claim that the algorithm will produce correct bounds - this we verify 
a posteriori in interval arithmetics. 


1. First we adjust C and s. Recall that, by (l42l) for k > M we want to choose C 
and s such that 


C _ aD 

¥ > 2 k s + l {j3~k- 2 ) ' 


(64) 


We want to increase s - therefore, we set the new parameters by 


s := s + 1, 


C := 


aD 

2(/3 — (M + l) -2 ) ' 


(65) 


2 . 


Trying to comply with ED we set the new u k s and u r k s 
1,... ,M by 



Nj + fc~ 2 /‘ 
2 {j3k 2 - 1) 
N r k + k- 2 r k 
2(/3fc 2 - 1) 


inductively for k = 


( 66 ) 


After each run we check inequalities (RTTl) and (1351) to see whether we obtained isola¬ 
tion. We may however require additional iterates to improve s. 


4.5. Proof of Theorem |l] 


For each of the parameter values and forcing terms we take a guess on the initial 
bounds with higher modes decay of order C /fc 4 . Once a guess for some given range of 
e is found it is easy to adapt it to another e by rescaling proportionally to the upper 
bound on |e|. After two iterates of procedure given in Subsection 14.41 we obtain the 
values of {u k }, C and s = 6 such that inequalities (1611) and (1551) hold. We present 
the approximate values (the first 5 significant digits of the actual values) in Table [3] 
We remark that in all of the cases it was enough to take M = 6. From Corollary (jHJ) 
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Tab 3.: Parameters M and s and approximate values of C, u l k and u r k , k = 1, ..., M 
used in the proof of Theorem Q] 


fe 

J T 

r 3 

J T 

P 

1.5 

1.75 

2.5 

1.5 

1.75 

2.5 

e 

[-0.05,0.05] 

[- 0 . 1 , 0 . 1 ] 

[-0.3, 0.3] 

[-0.05,0.05] 

[- 0 . 1 , 0 . 1 ] 

[-0.3,0.3] 

M 

6 

6 

6 

6 

6 

6 

u\ = u] 
«2 = ~ U 2 
“3 = ~ U 3 
U 4 = -U 4 

“5 = ~ U \ 
< = 

0.05743 
0.004018 
0.00022427 
1.1242 x 10 ~ 5 
5.7862 x 10 - 7 
5.5904 x 10 - 7 

0.082489 
0.0069984 
0.0004798 
2.9597 x 10 ~ 5 
1.9415 x 10 ~ 6 
1.9328 x 10 ~ 6 

0.16777 

0.020237 

0.0019934 

0.00017727 
1.8646 x 10 ~ 5 
2.1631 x 10 ~ 5 

0.059242 

0.0055667 
0.00054628 
9.3307 x 10 ~ 5 
2.9174 x 10 ~ 6 
7.7255 x 10 ~ 7 

0.085611 

0.0097091 
0.0010739 
0.000179 
7.6705 x 10 ~ 6 
2.596 x 10 ~ 6 

0.17515 

0.026475 

0.0035043 

0.00055121 
4.3434 x 10 ~ 5 
2.6004 x 10 ~ 5 

C 

4.6941 

13.039 

100.64 

4.8878 

13.613 

102.99 

S 

6 

6 

6 

6 

6 

6 


we conclude that for each Galerkin projection (1301) for n > M there exists a periodic 
solution (u + ’ n {t),u~' n {t)) such that 

M n 

« ±,n (t)e III"*’"*] ® n [~c/k',c/k a ], VteR. (67) 

k =1 k=M +1 

The above computations were done on a computer in interval arithmetics, as it 
would be tedious to do them by hand. The source files are available online [7' . The 
program uses interval arithmetics implementation from the CAPD package ill. 

By the change of variables (EHD we return now to the original coordinates Uk, Vk, 
i.e. the Fourier coefficients of u and Ut and obtain a sequence of periodic solutions 
(u n (t),v n (t)) of the Galerkin projections of the system ([25]) . From equation (1671) and 
the form of our change of variables it follows, that there exists a C > 0 (an exact 
value of which is not important to us), such that 

«£(*)<§, «*(*)<£■ (68) 

for all n, k : k < n and tet. A standard argument (cf. Theorem 10 in [2D]) proves 
that the set 

OO 

Y[ [- C/k 6 , C/k 6 } x [-C/k 4 , C/k 4 } (69) 

fc =i 

satisfies conditions C2, C3 and forms self-consistent bounds for (1251) . Note that 
at this moment we need the polynomial coefficient decay rate to be of order at 
least 2 for v k ’s and 6 for u k 's (we have 4 and 6, respectively). Let u n {t,x) = 
X^fceN u k(t) elkx , v n (t,x) = X^fceN v k(t) elkx . From Theorem [5] it follows that the 
sequence {(u n (t,x),v n (t,x))} n has a subsequence {(u ni (t, x), v ni (t, x))}i converging 
uniformly on compact time intervals to a solution (u*(t,x),v*(t,x)) of 


u t = v , 

'Ct = U X X 4” filLxxXX 4“ <j{u )xx 4” (^7 %) 7 

i.e. the Boussinesq equation m rewritten as a first order system. We have 
u*(t, x) = lim u ni (t, x) = lim u ni (t + r, x) = u*(t + r, x) 

l—too l—too 


(70) 


( 71 ) 
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for all hence the solution is periodic. 

The C° and L 2 bounds on u* and v* are computed from equations (l26l) and (l67l) . 
From the coefficient decay (lC>51) and elementary facts about the Fourier series (see 
Section 6 in 18|) it follows that u*{t,x) is of class C 4 and v*(t,x) is of class C 2 as 
functions of x. 
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